Mon. Not. R. Astron. Soc. 000, 1-?? (201 1) Printed 8 February 2012 (MN I*TgX style file v2.2) 



Systematic study of X-ray Cavities in the brightest galaxy of the 
Draco Constellation NGC 6338 



M. B. Pandge 1 , N. D. Vagshette 1 , L. P. David 2 , M. K. Patil 1 * 

1 School of Physical Sciences, Swami Ramanand Teerth Marathwada University Nanded, 431 606, India. 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden St. Cambridge, MA 02138 



8 February 2012 



ABSTRACT 

We present results based on the systematic analysis of currently available Chandra archive 
data on the brightest galaxy in the Draco constellation NGC 6338, in order to investigate 
the properties of the X-ray cavities. In the central ~6kpc, at least a two and possibly three, 
X-ray cavities are evident. All these cavities are roughly of ellipsoidal shapes and show a 
decrement in the surface brightness of several tens of percent. In addition to these cavities, 
a set of X-ray bright filaments are also noticed which are spatially coincident with the Ha 
filaments over an extent of 15 kpc. The Ha emission line filaments are perpendicular to the X- 
ray cavities. Spectroscopic analysis of the hot gas in the filaments and cavities reveal that the 
X-ray filaments are cooler than the gas contained in the cavities. The emission line ratios and 
the extended, asymmetric nature of the Ha emission line filaments seen in this system require 
a harder ionizing source than that produced by star formation and/or young, massive stars. 
Radio emission maps derived from the analysis of 1 .4 GHz VLA FIRST survey data failed to 
show any association of these X-ray cavities with radio jets, however, the cavities are filled by 
radio emission. The total power of the cavities is 17 x 10 42 erg s _1 and the ratio of the radio 
luminosity to cavity power is ~ 10 -4 , implying that most of the jet power is mechanical. 
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1 INTRODUCTION 

It is well established that feedback from active galac- 
tic nuclei (AGN) plays an important role in the evolution 
of hot gas in individual galaxies as well as galaxies in 
groups and clusters (Hoeft & Bruggen 2004, Dunn & Fabian 2006, 
McNamara & Nulsen 2007, Diehl & Statler 2008, David et al. 
2009, Gastaldello et al. 2009, Sun et al. 2009, Dong et al. 2010). 
Recent high-resolution X-ray observations with the Chandra and 
XMM-Newton telescopes have shown that the X-ray emitting gas 
from cooling-flow galaxies does not cool directly from the plasma 
phase to the molecular phase, but instead is reheated. The current 
paradigm is that AGN in the cluster dominant galaxy are respon- 
sible for this reheating, implying a complicated feedback process 
between the cooling gas and the central galaxy (Soker & Pizzolato 
2005). AGN also play important roles in shaping the morpholo- 
gies of the hot gas halos surrounding individual galaxies, groups 
and clusters and the formation of cavities and bubbles (Birzan et al. 
2004). These cavities are depressions in the X-ray surface bright- 
ness filled with low density relativistic plasma and are consequence 
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of the thermal gas being displaced by jets from the central AGN 
(Fabian et al. 2003). It is widely accepted that mass accretion onto 
the central massive black hole is the ultimate source of energy for 
the formation of such cavities, however,but there are still many de- 
tails left uncovered. 

Cavities embedded in hot X-ray halos are useful in quantify- 
ing the energy and power output of recent AGN outburst in central 
dominant galaxies. Present day observing facilities at X-ray wave- 
lengths with superb angular resolution have made it possible to 
detect such cavities or bubbles of a few to few tens of kpc size. 
Clusters of galaxies and have provided the strongest observational 
evidence for the AGN feedback in galaxy clusters. X-ray obser- 
vations with these telescopes and their systematic analysis have 
shown that about 20-25% of group and cluster galaxies harbour 
such cavities, while this detection rate may go even up to about 
70-75% in the case of X-ray bright cool-core clusters (Dunn et al. 
2005, Dong et al. 2010). In the majority of systems these X-ray 
cavities are found to be associated with radio jets and are spa- 
tially coincident with radio bubbles or lobes (McNamara & Nulsen 
2007). Combined X-ray and Radio observations of such systems 
have indicated that these cavities are inflated by bipolar jets ema- 
nating from AGNs located at the centre of a CDG (McNamara et al. 
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2009). Like radio jets, in many cases X-ray cavities also appear in 
pairs (Wise et al. 2007). These X-ray cavities are found to be filled 
with radio emission and in some cases are found to be connected 
with nuclei of host galaxies through synchrotron jets and tunnels in 
hot gas (Clarke et al. 2005). If these cavities are buoyantly driven 
outward, it is possible to estimate their dynamical age from the 
distance to the AGN. Though X-ray cavities have also been de- 
tected in groups, the majority of the studies of AGN feedback have 
been focused on cavity systems seen in the massive and X-ray lu- 
minous galaxy clusters, and groups have received less attention. In 
fact, most of the baryonic matter in the Universe resides in groups 
and hence they are perfect laboratories to understand the impact of 
feedback on they formation and evolution of galaxies. 

The centrally dominant galaxies (CDGs) that lie at centres 
of cooling-flows are often found to harbour cool interstellar me- 
dia which can be traced by nebular emission lines (Heckman et al. 
1989, Crawford & Fabian 1993), neutral hydrogen (O'Dea et al. 
1994,Taylor 1996), molecular gas (Edge et al. 2002), and associ- 
ated star formation (Crawford et al. 1999,Rafferty et al. 2008). Far- 
UV observations combined with optical and IR data can constrain 
the star formation history and the properties of young stars associ- 
ated with emission line nebulae in these cooling flow galaxies. Evi- 
dence for star formation in cool core cluster CDGs have been gath- 
ered from optical and UV observations (Crawford & Fabian 1993, 
Wise et al. 2004, Rafferty et al. 2006). A strong correlation has 
been noticed between the occurrence and strength of star forma- 
tion rates in CDGs and their X-ray cooling rates measured with the 
Chandra and XMM-Newton observatories (McNamara & Nulsen 
2007). However, in contrast to the cooling rates of several hundreds 
of M©yr -1 in cooling flow clusters , star formation rates are only 
about a few to a few tens of M yr _1 . 

In this paper, we present a detailed analysis of the X-ray data 
on NGC 6338, a brightest galaxy in the Draco constellation, to in- 
vestigate the physical properties of the X-ray cavities hosted in this 
system. NGC 6338 is the brightest member in the compact group 
of more than 13 galaxies within a field of 8 x 8arcmin 2 , and forms 
a close pair with MCG+10-24-117 at l'.2 N. Other members in- 
clude NGC 6345 at 3'.6 S, NGC 6346 at 5'.3 S and IC 1252 at 
about 4'. 5 SE. NGC 6338 is considerably more luminous than the 
other members having an absolute magnitudes of Mb = —22.2, 
M v = -21.92 (Martel et al. 2004) and is well described by the 
r 1//4 surface brightness profile law. The NGC 6338 group seems to 
be relaxed, with a smooth velocity distribution about a mean red- 
shift of cz = 8222 km s _1 (Tabel.l). The NGC 6338 hosts a central 
bright source as seen in the X-ray emission. NGC 6338 is a previ- 
ously studied object and offers one of the best chances of studying 
interaction between the active galactic nucleus (AGN) and inter- 
stellar medium (ISM) in its neighbourhood. HST observations of 
this system have detected asymmetric Ha filaments extended up to 
about 7.5 kpc (Martel et al. 2004) and dust in central regions. 

The structure of this paper is as follows: Section 2 describes 
the Chandra observations and data preparation. Section 3 discusses 
the imaging analysis and the detection of the X-ray cavities as well 
as their spectral analysis. Section 4 presents a comparison of the 
X-ray data with the optical and radio data. Calculations of the en- 
ergetics associated with the central AGN and the results are sum- 
marised in Section 5. All distance dependent parameters have been 
computed assuming the luminosity distance of NGC 6338 equal to 
1 15.3 Mpc and H =70 km s' 1 Mpc -1 . 
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Figure 1. Raw 0.3-3.0 keV ACIS-I Chandra X-ray image of the central 
5xAarcmin 2 region of the brightest group galaxy NGC6338 overlaid on 
the optical D25 ellipse. The arrows indicate the features like, cavities, fila- 
ments, etc. as discussed in the text (Section 3). A nearby bright member at 
about 42.7 kpc is also shown in this figure. 



Table 1. Global parameters of NGC 6338 



RA& DEC 
Morph 

Mag(B T ) 

M B 

Size 

Distance(Mpc) 
Redshift(z) 

Radial Velocity (km s" 



17:15:23.0; +57:24:40 

SO 

13.6 

-22.2 

l'^xl'.O 
115.3 
0.02742 
8222 



2 OBSERVATIONS AND DATA PREPARATION 

NGC 6338 was observed by the Chandra X-ray Observatory on 
September 17-18, 2003 for an effective exposure time of 47.94 ks 
(ObsID 4194). NGC 6338 was located on the 13 chip of the ACIS-I 
detector. We have reprocessed X-ray observations using the stan- 
dard tasks available within CIAO 1 4.1.0 and CALDB 4.1.2 pro- 
vided by the Chandra X-ray Centre (CXC), and following stan- 
dard Chandra data-reduction threads 2 . Periods of high background 
flares were identified using 3 a clipping of the full-chip light using 
the Icjclean task and binned in the lengths of 200 s. After remov- 
ing these periods of high background, the cleaned data have a net 
exposure time of 44.5 ks. Point sources on the 13 chip were identi- 
fied using the CIAO task wavdetect with the detection threshold of 
10~ 6 . 



3 RESULTS 

3.1 Investigation of X-ray cavities in NGC 6338 

A raw 0.3-3.0 keV ACIS-I Chandra image of NGC 6338 is shown 
in Figure 1. The nearby member MCG +10-24-117, separated by 
42.7 kpc from NGC 6338, is also shown Figure 1. From this fig- 
ure, it is clear that the hot gas in the central region of this galaxy 
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Figure 2. Left panel: 0.3-3.0 keV ACIS-I Exposure corrected X-ray image of the central region of NGC 6338. Middle panel: Unsharp-masked image of NGC 
6338 in the energy band 0.3-3.0 keV produced by subtracting a 5a wide Gaussian kernel smoothed image from that smoothed with 2a Gaussian kernel. Right 
panel: 0.3-3.0 keV elliptical 2D beta model subtracted residual image. Both these figures clearly reveal at least a pair of cavities (shadow regions marked as 
NE and SW cavities) and filaments in the form of excess emission (white shades). 




Figure 3. Left panel: Azimuthally averaged, clean background subtracted 0.3-3.0 keV surface brightness profile of NGC 6338. The best fit (3 model is shown 
as a solid line. Right panel: Clean background subtracted 0.3-3.0 keV surface brightness profiles of X-ray photons extracted from conical sectors covering NE 
(red colour) and SW (blue colour) cavities. For the NE cavity photons were extracted between 110° — 170°, while that for SW cavity photons were extracted 
from 270° — 340° (both measured from west to north). Depressions in the profiles are visible at the location of cavities. 



is not distributed smoothly but instead shows fluctuations and dis- 
continuities. A pair of cavities can be seen in this Figure and are 
identical to those reported by (Dong et al. 2010). However, to visu- 
alize and study these features in detail we have derived an unsharp 
masked image as well as a residual image after subtracting a 2D el- 
liptical f3 -model (Dong et al. 2010). The methods methods behind 
the derivations are discribed in more detail below. 

3.1.1 Unsharp -Masked Image 

To enhance the visualisation effect of the faint surface brightness 
fluctuations seen in the NGC 6338 (see Fig. 2 Left-hand panel) 
at fine spatial details, we have produced a 0.3-3.0 keV unsharp- 



masked image. For this purpose the exposure-corrected, back- 
ground subtracted image was first smoothed with a narrow Gaus- 
sian kernel of width 2a (\a=\ pixel) using the CIAO task acon- 
volve, so that it suppresses pixel-to-pixel variations while preserv- 
ing structures on likely scales. A similar image was also generated 
by smoothing the data with a wider Gaussian kernel of 5 a, so that 
it preserves overall morphology of the host galaxy while erasing 
small-scale features. The unsharp-masked image was the generated 
by subtracting the image smoothed with the 5a from that smoothed 
with the 2a Gaussian (see Dong et al. (2010) for details). The re- 
sulting unsharp-masked image is shown in Figure 2 (Middle panel) 
which reveals at least two X-ray cavities, one on the north-east (NE) 
side and other on the south-west (SW) side of the optical centre of 
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NGC 6338. Both of these cavities are roughly of ellipsoidal shapes 
and are located at a projected distance of about 3.2 kpc from the 
centre of the host galaxy. Apart from these two prominent cavities, 
several other surface brightness discontinuities or "holes" and "fil- 
aments" are also evident in this figure. One more relatively fainter 
cavity-like structure is also visible in this figure at a distance of 
about 10 kpc in the northern (N) direction from NGC 6338. 



3.1.2 Elliptical 2D f3 model 

The features seen in the unsharp-masked image were highlighted 
using a residual map of NGC 6338 after subtracting its smooth 2D 
elliptical beta model from the background subtracted image. For 
this purpose, an elliptical 2D f3 model was fit to the clean back- 
ground subtracted image of NGC 6338 in the energy range 0.3- 
3 keV using Sherpa available within CI AO. The fitting parameters 
i.e., ellipticity, position angle, normalisation factor and local back- 
ground values, were kept free during this fit. The resulting best-fit 
model with f3 = 0.89 ± 0.02 and r c =9.10±0.40 kpc was then sub- 
tracted from the exposure corrected image of NGC 6338 to produce 
a residual image Figure 2 (Right panel). This image reveals several 
departures from the averaged X-ray emission due to asymmetry in 
hot gas morphology. These departures appear either in the form of 
depressions at the locations of cavities (SW and NE cavities) or ex- 
cess emission in the form of bright filaments (along NW and SE). 
Apart from two prominent cavities, a pair of holes or bubbles are 
also evident in this figure near the NW and SE filaments. 



3.1.3 Surface Brightness Profile 

A surface brightness profile of the extended emission from NGC 
6338 in the energy range 0.3-3.0 keV was derived by extracting 
counts from concentric annuli up to a radius of 30 // from NGC 6338 
and is shown in Figure 3. From figure, it is evident that the X-ray 
emission from NGC 6338, like in other cooling flow group galax- 
ies, shows a central peak in the X-ray surface brightness profile. 
Depressions in the surface brightness profile, due to presence of 
cavities and a small excess emission at about 10" due to filaments, 
are evident in figure. Despite the fluctuations due to cavities and 
filaments, the azimuthally averaged surface brightness profile of 
NGC 6338 is remarkably smooth and the standard one-dimensional 
P model can be fit to this profile 
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where ro is the core radius. Best-fit parameters are ro = 
2.99 =b 0.34 kpc and f3 = 0.46 =b 0.014 and are in good agreement 
with those reported by Dong et al. (2010). A careful look at the 
azimuthally averaged surface brightness profile reveals a marginal 
break at about 30 kpc from the centre of NGC 6338 which may be 
a cold front due to sloshing Markevitch & Vikhlinin (2007). 

To highlight, the depressions seen in the X-ray emission from 
the cavities, we have generated a 0.3-3.0 keV surface brightness 
profile of the background subtracted, X-ray image of NGC 6338 by 
extracting photons from two conical sectors, Sector 1 (covering NE 
cavity) and Sector 2 (covering SW cavity). The angular coverage of 
these sectors in west to north directions are 110° to 170° and 270° 
to 340°, respectively. The resultant surface brightness profiles are 
shown in Figure 3 (Right panel) and show noticeable azimuthal 
variations. Apart from two main cavities, several other azimuthal 
variations in the surface brightness between about 5"and 30"of are 




Figure 5. Left Panel: 0.3-3.0 keV background subtracted clean X-ray image 
of NGC 6338 along with an annular ring, covering both cavities and fila- 
ments, divided into 14 different sectors. Right Panel: Profile showing vari- 
ations in count rates extracted from these 14 sectors. Sectors corresponding 
to cavities show significantly less counts, while those corresponding to fil- 
aments show excess counts. The horizontal dashed line represent expected 
counts in the absence of such deviations. 



evident and may indicate sloshing motions of hot gas in NGC 6338. 
To investigate these features in more detail, we have derived surface 
brightness profiles showing variations in two different strips shown 
in Figure 4, where strip 1 covers both cavities while strip 2 covers 
filament regions. At the location of the NE cavity a 10% deficiency 
is seen along with a 20% deficiency at the location of SW cavity. 
An overall excess was seen in the strips covering the filaments. 

Figure 5 (Left panel) shows an annular ring, covering both the 
cavities as well as the filaments, divided into 14 different sectors. 
0.3-3.0 keV X-ray counts extracted from each sector were plotted 
as a function of the sector number and are shown in Figure 5 (Right 
panel). This profile shows two peaks between sector numbers 1- 
4 and 8-11, while two depressions are also seen between sector 
numbers 5-8 and 11-14, respectively. Peaks and depressions in this 
figure coincide with the filaments and cavities, seen in the unsharp- 
masked image of NGC 6338. The horizontal dashed line in this fig- 
ure represents the expected average counts in the absence of such 
deviations. A small depression is also noticed between sector num- 
bers 3-5 which represent the partially covered third cavity in the 
north (N) direction. 



3.2 Spectral Analysis 

The properties of the complex X-ray morphology due to bubbles 
and cavities in NGC 6338 cannot be fully investigated from sim- 
ple, azimuthally averaged surface brightness profiles. Therefore, to 
derive the physical properties of these features, we have performed 
a spectral analysis of X-ray photons extracted from these different 
regions. 



3.2.1 Projected gas temperature profile 

To establish the general radial dependence of gas properties, we 
have extracted spectra from the 13 chip in concentric circular an- 
num' in the 0.5-5.0 keV energy band. These annulii were centered 
on the peak of X-ray emission, which corresponds to the optical 
centre of NGC 6338. Source spectra, background spectra, photon 
weighted response files, and photon weighted effective area files 
were generated for each region using the CIAO task acisspec. The 
spectra extracted from each annulus were initially fitted with a sin- 
gle temperature MEKAL model. However, better results were ob- 
tained when the spectra were fit with an absorbed two tempera- 
ture vapec thermal plasma model with the neutral hydrogen col- 
umn density fixed at the Galactic value, Nh = 2.60 xlO 20 cm -2 
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Figure 4. Left panel: 0.3-3.0 keV unsharp masked image showing two rectangular strips. Strip 1 covers both cavities while strip 2 covers filament regions. 
X-ray photons from these strips were extracted to visualise surface brightness fluctuations in a better way. Middle panel and Right panel: Profiles showing 
variations in count rates extracted from strip 1 and strip 2, respectively. Locations corresponding to cavities show significant deficit of counts whereas locations 
corresponding to filaments show excess of counts. Peaks in these profiles corresponds to centre of NGC 6338. 
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Figure 6. Unsharp masked image showing different regions from where 
X-ray photons were extracted for studying their spectral properties. 



Rasmussen & Ponman (2007). This fitting allowed us to constrain 
the abundances of O, Ne, Mg, Si, S, and Fe simultaneously, in addi- 
tion to the temperature and normalisation. The second vapec com- 
ponent was added to account for the emission originating from gas 
at multiple temperatures. We also added a power-law component 
of photon index fixed at r=1.4 to account for emission from unre- 
solved LMXBs. 

The azimuthally averaged X-ray temperature profile is shown 
in Figure 7 (Left panel). From this profile, it is clear that 
the temperature of the X-ray emitting gas increases monotoni- 
cally from 1.51±0.14keV at centre to 3.39±1.45keV at about 
50" which is similar to that seen in other cooling flow galaxies 
(Gastaldello et al. 2007, Rasmussen & Ponman 2007, David et al. 
2009, Sun et al. 2009). We have also derived electron density and 
pressure profiles (Figure 7 middle and right panels ), respectively. 
Rasmussen & Ponman (2007) derived a radial temperature profile 
for NGC 6338 by fitting single-temperature models over the radial 
distance of about 657 kpc and have shown that the temperature of 
hot gas in this system increases monotonically from about 1.29 keV 
at the core region to a maximum of 3. 1 keV at 65.7 kpc and then de- 



creases outwardly. Our estimates of the temperatures with in 26 kpc 
are in good agreement with those derived by Rasmussen & Ponman 

(2007) . They also defined this system as the hottest system in their 
sample and found a metal excess well beyond the central region. 
Our estimates of other thermodynamical parameters from this anal- 
ysis are also found a match well with the estimates of Rafferty et al. 

(2008) . 



3.2.2 Cavities and Filaments 

To deduce the spectral properties of the hot gas in the cavities and 
filaments in more detail, we extracted spectra representative of both 
the cavities (NE and SW) as well as of both filaments (NW and SE) 
using elliptical apertures (yellow and blue colour ellipses shown in 
Figure 6). An absorbed apec model was fit to photons in the 0.5- 
5.0 keV energy band and yield temperatures of 1.72±0.11 keV and 
1.72±0.08 keV, for the NE and SW cavities, respectively. The spec- 
tra extracted from the NW and SE filaments were also fit with an 
absorbed apec model, which yield temperatures of 1.28±0.04keV 
and 1.58±0.09keV for NW and SE filaments, respectively. This 
means gas associated with the cavity regions is hotter than that from 
the filament regions and is perhaps due to the emission from fore- 
ground and background gas. We have estimated temperature of the 
undisturbed gas at the same radius as the filaments by fitting ab- 
sorbed apec model to the spectrum extracted from this region and 
was found to be equal to 2.21 ±0.20 keV, which is again higher than 
the filament temperatures. 

To investigate the nature of these brighter filaments with softer 
X-ray emission, we plot the spectra extracted from both the fil- 
aments as well as the cavities in Figure 8. From this figure it is 
clear that spectrum taken from the NW filament has a more promi- 
nent Fe-L feature which is shifted towards lower excitation energy 
(green data points) compared to those from the cavity regions (red 
and black data points). This confirms that the gas in the filaments 
is cooler than the gas in the cavities. These cooler filaments are co- 
inciding with the Ha emission lines filament and dust knots (Fig- 
ure 10). 
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Figure 7. Azimuthally averaged projected profiles of temperature (left), electron density (middle), and pressure (right) measured with Chandra. 



Energy (keV) 

Figure 8. Spectra extracted from NE cavity (black data points), SW cavity 
(red data points) and NW and SE filaments (green data points). Overall 
shape of the Fe-L hump in the case of filaments is significantly different and 
is shifted toward lower temperature side, indicating that hot gas in filaments 
is relatively cooler than that contained in cavities. 



3.2.3 The Central Source 

A 2-7 keV X-ray image of NGC 6338 reveals a bright X-ray point 
source located at its centre, which coincides with the radio source 
in the FIRST image and also with the optical centre of the galaxy. 
To examine the spectral properties of probable AGN, we extracted 
a 0.5-7 keV spectrum from within the central 2 // (^1.18kpc). We 
tried to fit this spectrum with a power-law component, however, it 
resulted in poor fit. Therefore, a combined thermal and power-law 
component was used. The best-fit temperature in this double com- 
ponent model is 1.49 keV and is consistent with the central gas tem- 
perature for this galaxy. The best-fit power-law index in the present 
case was found to be T=2.6, which signifies that the central source 
has a sufficiently hard component like typically seen in AGN can- 
didates. The hydrogen column density was kept as a free parameter 
during this fit, however, we do not find significant excess of it rela- 
tive to the Galactic value. From this thermal plus power-law model, 
we estimate the unabsorbed 0.5-7 keV flux of the central source to 
be 6.52xl0 -14 erg s _1 cm -2 leading to a luminosity of 1.12x 
10 41 ergs- 1 .(Tabel.2 




Figure 9. 0.5-8.0 keV X-ray spectrum extracted from the 30"region cen- 
tered on NGC 6338. The continuous line represents the best-fit model. Cen- 
tral 2 // region is excluded in this analysis. 



3.2.4 Total Diffuse Emission 

In order to determine the average properties of the diffuse gas NGC 
6338, we extracted a spectrum from within 30 // (16.2kpc), exclud- 
ing the central 2 // region as well as other point sources. A total of 
11855 background- subtracted counts were extracted in the range 
between 0.5-7 keV. We initially fitted the spectrum with a single 
apec model, however, the resulting model showed residuals near 
Si (1.82 keV), S (2.38 keV), Ar. (3.06 keV) and Fe (6.3 keV) lines. 
Therefore, we use the vapec model by allowing abundances of O, 
Ne, Mg, Si, S, and Fe to vary independently. A much better fit was 
obtained with abundances Si (~ 0.61 ± 0.14), S (~ 0.55 ± 0.17), 
Ar (~ 0.37 ± 0.37), and Fe (~ 0.49 ± 0.09) and an average tem- 
perature of 1.49±0.17keV. To account for the contribution from 
unresolved point-like sources, we also included a power-law com- 
ponent with r=1.4. (see Fig.9 and Tabel.2) 



4 DISCUSSION 

4.1 Optical and UV analogue of the X-ray emission 

A systematic study of cooling flow clusters using high resolu- 
tion data from the Chandra and its optical counter part from 
the HST have confirmed the strong spatial correspondence be- 
tween X-ray filaments and optical emission line (Ha ) filaments 
(McNamara et al. 2000, Fabian et al. 2001, Blanton et al. 2001, 
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Figure 10. From left to right: (1) HST Ha: emission line filaments (white shades) noticed in NGC6338. (2) Ha emission line filament contours (green 
contours) superposed on the 0.3-3.0 keV X-ray unsharp mask image of NGC 6338. Note the spatial correspondence of these filaments with the filaments in 
X-ray emission. Dust contours (red contours at the centre) are also displayed for comparison. (3) HST I-band (F814W) optical image of NGC 6338. Embedded 
dust features are evident in this figure. (4) 2D elliptical model subtracted residual I band image of NGC 6338 (magnified image). Dust knots (white shades) 
are clearly visible in the central part of this galaxy. 



Table 2. Spectral analysis of cavities, filaments, total diffuse component 
and the central source 



Regions 


N H 
(10 20 cm" 2 ) 


Best-fit 


kT 
(KeV) 


Goodness of fit 
X 2 /dof 


NE Cavity 


(2.14) 


apec 


1.72±0.11 


1.14 


SW Cavity 


(2.14) 




1.72±0.08 


0.75 


NW Filament 


(2.14) 




1.28 ±0.04 


0.87 


SE Filament 


(2.14) 




1.58±0.09 


1.16 




(2.14) 




2.03 ±0.42 


0.75 


Central source 


(2.14) 


apec+powlaw 


1.01 ±0.09 


1.24 


nponent 


(2.14) 


vapec+powlaw 


1.50±0.05 


1.14 



Notes on column: Col.l- Regions used to extract the spectra in the energy band (0.3-5.0 keV), col. 2- 
Galactic hydrogen column density fixed during spectral fitting, col. 3- Best-fit model, col.4-Temperature of the hot gas 
col.5-Goodness of fit 



Young et al. 2002). To examine a similar association between the 
two filaments, we have made use of the HST archival data using 
the filters F814W and FR656N with integration times of 700 s and 
1700 s, respectively. The details of the observations and data reduc- 
tion can be found in Martel et al. (2004). From the analysis of these 
observations, it is revealed that, the cooling core of NGC 6338 har- 
bours a remarkable Ha +[N/i] emission line filamentary system 
which is oriented along the NW to SE direction with a total exten- 
sion of about 15 kpc (Figure 10-1) and exhibits a surprising corre- 
spondence with the X-rays image. Unlike other systems, both the 
Ha and X-ray filaments are oriented perpendicular to the X-ray 
cavities seen in this galaxy. Further, more the Ha emission line fil- 
aments in this galaxy are not continuous but are broken. Filaments 
oriented in the NW direction are relatively fainter and diffuse, while 
those in the SE direction are divided into two arms, one oriented 
along the south and other along the east. Contours of Ha filaments 
are overlaid on the X-ray emission map (Figure 10-2) which reveals 
a strong spatial correspondence with X-ray filaments. 

In cooling flow systems, the optical emission line filaments 
are often found to be associated with the interstellar dust. Interstel- 
lar dust is also evident in between the central region of the HST I 
band image of NGC 6338 (Figure 10-3). To assess the association 
between the dust and ionized gas in NGC 6338, we have derived 
residual map of NGC 6338 after subtracting its best fit 2D ellipti- 
cal de Vaucouleurs model (Patil et al. 2007) from the I band image 
(Figure 10-4). The residual map derived for this system exhibits 
the presence of a disturbed dust morphology in the form of fila- 
ments or knots (white shades) in the central region of this galaxy 
and is found to be extended up to about 3.6 kpc. These results are 
in full agreement with those reported by Martel et al. (2004). The 
dust knots seen in the southward direction correspond with the Ha 



emission line filaments as well as X-ray filaments in the central re- 
gion, however, the dust extent is very small compared to the Ha and 
X-ray emission. The brightest dust knot in the northward direction 
of the extinction map coincides with the third, relatively fainter, de- 
pression seen in the X-ray emission map, implying that this may be 
due to the absorption from the dust patch. Thus, highly asymmetric 
filaments of Ha emitting gas have been found within the central 
~8.94kpc of NGC 6338 and are found to be associating with the 
optically absorbing dusty gas as well as soft excess emission fila- 
ments seen at X-ray wavelengths and provides a strong evidence 
for the hypothesis that dust has been buoyantly transported from 
the galactic cores out to several kilo-parsecs following a feedback 
heating events (Gastaldello et al. 2009). The extended morphology 
of Ha emission seen in this galaxy and its spatial correspondence 
with the X-ray filaments suggest that the optical emission probably 
does not arise from shocked condensation (see David et al. (1988)). 
We tried to check its correspondence with the Near Ultra Violet 
(NUV) data from the GALEX observatory, however, due to poor 
signal in the NUV image we could not assess such an association. 

In the recent imaging survey of brightest cluster galaxies with 
the Spitzer Space Telescope, Quillen et al. (2008) have defined this 
galaxy as an excess IR emission galaxy. The observed extended 
emission at mid-IR frequencies by Spitzer support the hypothesis 
that the dust has been buoyantly transported from the galactic core 
out to several kpc due to the AGN feedback heating mechanism. 
The disturbed morphology of the dust within the central 3.6 kpc 
region seen in the HST image provides additional evidence regard- 
ing a recent energy release from the central source. The impact of 
this energy release has already been evidenced through high resolu- 
tion X-ray imaging observations with the Chandra telescope in the 
form of a pair of cavities on either sides of the nucleus. A system- 
atic study of these cavities allows us to shed light on the amount of 
energy deposited by the central AGN and the time of such a release. 

4.2 Star formation 

In the absence of central heating, gas in the ICM cools below the 
X-ray temperatures and is gets deposited onto the central galaxy. 
The cooled ICM accumulates at the centre of the potential well in 
the form of molecular clouds which may lead to formation of stars 
. This has been confirmed through the common occurrence of gi- 
ant Ha filaments and on-going star formation in cool-core clus- 
ter BCGs (Crawford et al. 1999, Rafferty et al. 2008). Ha emis- 
sion maps derived from the analysis of HST archival data on NGC 
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Figure 11. Optical SDSS spectrum of NGC6338. Regions of O II emission 
lines and Ha; +[N II] lines are highlighted in lower panel. 

6338 galaxy have exhibited giant, asymmetric Ha emission line 
filaments extending to large distance with a spatial correspondence 
with the X-ray filaments (Figure 10). Spectroscopic analysis of the 
hot gas extracted from these filamentary structures reveals that the 
gas corresponding to the filaments is cooler than that in its vicinity 
as well as in cavities. This means that the cooler gas in the filamen- 
tary structures is cooling which may result in to the formation of 
stars. 

We have made an attempt to assess any ongoing star formation 
in NGC 6338 using optical spectroscopic observations available in 
the archive of the Sloan Digital Sky Survey 3 (SDSS) and is shown 
in Figure 11. Using the de blending technique of IRAF 4 , we have 
measured flux of the Ha emission lines from this galaxy. This en- 
abled us to measure the total Ha luminosity and estimate present 
star formation rate using relation given by Kennicutt (1998) assum- 
ing that all of the Ha emission arises from star formation. 

SFR = 7.9 x 10" 42 x LhcMq yr' 1 (2) 

where Lna is in units of erg s _1 . From the measured value of 
the Ha luminosity for this target galaxy we estimate the star for- 
mation rate (SFR) to be equal to 0.05 M© yr~ x , which is smaller 
than the estimates of 0.13M©yr _1 by Crawford et al. (1999). The 
discrepancy in the two estimates is because we have not corrected 
the Ha flux for intrinsic extinction. Crawford et al. (1999) gives an 
Ha emission line luminosity of 1.2 x 10 40 erg s _1 and defined this 
system as a low luminosity emission-line system. They estimated 
the line ratios [NII]/Ha - 3.6, [SII]/Ha - 0.46, [0III]/H/3 ~ 
0.78, from which they demonstrated that the source required for 
maintaining such a high level of ionisation must be hard and should 
occupy a smaller region of the host galaxy. To examine the possibil- 
ity of such an ionization by radiation field from young, high-mass 
stellar population, Fabian et al. (2011) estimated the population of 
05 stars to be equal to 1.2 x 10 5 and noticed that this number is not 
enough to achieve the level of observed emission line flux. Other 
possible sources for such ionisation mechanism can be either turbu- 
lent mixing of layers or a low-level nuclear activity. In the case of 

3 SDSS is the Data base http://www.sdss.org 

4 Image Reduction And Analysis Facility,The IRAF is a distributed by the 
National Optical Astronomy observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomyjnc, Under cooperative 
agreement with the National Science Foundation 




Figure 12. The radio contours from the 1 .4 GHz VLA FIRST observation 
(green colour) are overlay on the unsharp masked X-ray image of NGC 
6338. The radio contours are at 0.3, 0.8, 2.12, 5.65, 15.1 and 40 mJy/beam 
levels, image rms is ~ 140 /i Jy beam- 1 and beam FWHM is 5"(black 
circle). 

galaxies where harder ionisation are apparent in smaller regions, 
low level nuclear activity like Low -Ionisation Nuclear Emission 
line Regions (LINER), can be the appropriate source of ionization. 

4.3 Association with the radio source 

The X-ray cavities found in the intra-cluster medium (ICM) are 
thought to be inflated by ambipolar jets originating from the cen- 
tral AGN that injects energy into small regions at their terminal 
points and expand until they reach pressure equilibrium with the 
surrounding ICM (Blandford & Rees 1974). This results in the for- 
mation of a pair of under dense bubbles on the opposite sides of 
the nucleus. The majority of central dominant galaxies (CDGs) are 
believed to host a central radio source and hence an active galactic 
nucleus. 

To find a similar behaviour at the centre of NGC 6338, we 
looked for the availability of radio data in the NVSS and FIRST 
catalogues and found a detection at 1.4 GHz in the FIRST Faint 
Images of the Radio Sky at Twenty -Centimetres survey. Using this 
data set, we derived a radio emission map for NGC 6338, contours 
of which are overlaid on the unsharp-masked X-ray image of NGC 
6338 (Figure 12). This figure confirms the presence of a point-like 
radio source at its center with radio contours extending up to 6 kpc 
covering both cavities, and are consistent with those reported by 
Dong et al. (2010) and Jetha et al. (2007)). However, deep, high- 
resolution, multi frequency radio observations from GMRT are re- 
quired for assessing extended nature of this source and will be pre- 
sented in the future paper. 

4.4 Cavity Energetics 

Chandra images of groups and clusters hosting central radio 
sources have amply demonstrated that the X-ray emitting gas in 
the vicinity of AGNs is displaced due to interactions between hot 
plasma and the central source and results in surface brightness de- 
pressions like those seen in the X-ray emission maps of NGC 6338. 
These X-ray surface brightness depressions are cavities or bubbles, 
devoid of gas at the local ambient temperature (McNamara et al. 
2000), and can be used as a calorimeter for the nuclear activity. The 
kinetic energy of jets emanating from a central black hole can be 
quantified by knowing how much energy is required to inflate the 
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Table 3. Cavity Properties of NGC6338 



Cavity Parameters 


NE cavity 


SW cavity 


a (kpc) 


6.74 


6.20 


b (kpc) 


1.73 


2.80 


R (kpc) 


3.56 


4.10 


tc s (yr) 


0.57 xlO 7 


0.61 xlO 7 


tbuoyiy?) 


0.86xl0 7 


1.09X10 7 


trefilliy?) 


l.OOxlO 7 


1.26xl0 7 


P cavity (ergs -1 ) 


6.76 xlO 42 


10.25 xlO 42 



X-ray cavities (Birzan et al. (2004), Rafferty et al. (2006)). This es- 
timate along with age of the cavity can be used to quantify the 
power of the central engine. 

From the analysis of Chandra observations of this galaxy we 
have detected two, probably three, X-ray cavities in the ICM asso- 
ciated with this central galaxy. These cavities are apparently con- 
nected to the central source and are perhaps undergoing expansion. 
To estimate the kinetic power of the central engine responsible for 
the formation and inflation of these cavities, it is required to quan- 
tify the energy involved in creating such bubbles. For the case of 
slowly rising bubbles it is equal to its enthalpy and is sum of the 
internal energy within cavity and pdV work done by the AGN on 
X-ray gas 

E bubble = —^—rpV +pV= -^- P V (3) 
7 — 1 7 — 1 

where p is gas pressure and is estimated from analysis of X- 
ray data, V is volume of cavities and 7 is mean adiabatic index of 
the fluid within the bubble (Birzan et al. 2004, Rafferty et al. 2006, 
Dunn & Fabian 2006, Dunn & Fabian 2006, Allen et al. 2006, 
David et al. 2009). Assuming that the fluid in the cavities is rela- 
tivistic plasma (7 = 4/3), the total energy contained in the bubble 
is given by E 'bubble = 4pV '. As these estimates are based on mea- 
surements of the X-ray cavity sizes and surrounding gas pressures 
they are independent of the radio properties, and hence are the more 
reliable estimates. To estimate the power of the X-ray cavities ob- 
served in the cooling flow system NGC 6338, we have followed 
the procedure outlined by Birzan et al. (2004). In estimating the 
volume of these cavities, we assume they are prolate ellipsoidal 
shapes, with semi-major and semi-minor axes equal to a and 6, re- 
spectively, and V = 47ra 2 b/3. This leads to volumes of 3.18 x 10 65 
and 2.42 x 10 66 cra 3 , for NE and SW cavities respectively. 

To quantify the rate at which the energy is injected into the 
ICM by the central AGN, it is required to know the ages of cav- 
ities. Here the age of each cavity was estimated in three different 
ways, (i) by calculating the time required for the cavity to rise to its 
present location from the radio core at the speed of sound, (ii) by 
calculating the time required for the cavity to rise buoyantly at its 
terminal velocity, and (Hi) by calculating the time required to refill 
the displaced volume of bubble when its rises upward (Birzan et al. 
2004). The time required for a cavity to rise to its present distance 
from the central AGN with the speed of sound is given by, 

t Cs =- = Ry/nmH/l kT (4) 

c s 

where 7 = 5/3, /i = 0.62 and R is the projected distance from 
centre of the cavity to the central AGN. With time the cavity may, 
detach from the AGN and hence rise buoyantly through the ICM to 
its present location on a time scale (Churazov et al. 2001) , 

huoy = R^/SC D /2gV (5) 
where Cd = 0.75 is the drag coefficient, S (= tto 2 ) is the cross- 



sectional area of the cavity and g = GM/R 2 is the gravitational 
acceleration (Birzan et al. 2004). Finally, the time required by the 
material to refill the displaced volume as the bubble rises upward is 
given by, 

Uefui = 2R^r/GM = 2^f7fg (6) 

where r is radius of the cavity and is equal to \fab for an 
ellipsoidal shape, In the present study we have estimated the ages 
of both cavities using all three methods discussed above and are 
given in Table 3. From the comparison of the three estimates it is 
clear that t Cs <C tbuoy <C t re fiii and all are within a factor of two. 

Here t Cs represents the time required for the cavity to expand 
to its present volume and is appropriate for active cavities, where 
cavities are powered by the central AGN, while tbuoy measures the 
rise time of cavities and is appropriate in the case where cavities 
are detached from the AGN and are buoyantly rising through the 
ICM i.e., the cavities which are not associated with the radio jets 
(ghost cavities). In the present case of NGC 6338, both cavities are 
well detached from the AGN and are filled in by radio emission at 
1.4 GHz. However, our radio analysis failed to detect any signatures 
of association of the cavities with the radio jets. Therefore, it is not 
clear which estimate suits well the age of the cavities in NGC 6338. 
As a result, the average of three estimates was taken as the age of 
cavities and was used for quantifying their power using 

Pcavity — -E'totai/ ^ tage ^ (7) 

and is to 17.01±3.25 x 10 42 erg s _1 . If we add the power due to the 
third cavity then this estimate could be even larger. As these esti- 
mates do not include possible hydrodynamical shocks Pcavity rep- 
resents a lower limit to the total power of the AGN. Cavagnolo et al. 
(2010) have also quantified the total power of both cavities through 
a similar analysis and have found it to be ll.O^g x 10 42 erg s _1 . 
Uncertainties in the two estimates arise because they used buoy- 
ancy rise time as age of the cavity. Though, all of the time esti- 
mates in the present study are within a factor of two, it leads to the 
reasonable uncertainty in the cavity power. Moreover, uncertainty 
in the volume of the cavities will amplify this discrepancy, even if 
we neglect uncertainties in the temperature, pressure and electron 
density. 

A comparison of total power with the observed X-ray lumi- 
nosity within the region occupied by cavities shows that total cav- 
ity power can balance less than one-half of radiative losses within 
the central 10 kpc as in NGC 5044 (David et al. 2009). The cool- 
ing time for NGC 6338 cluster is - 0.5 x 10 8 yrs (Rafferty et al. 
2008) and is sufficiently larger than the estimated age of cavities. 
High resolution and high sensitivity radio observations at 1.4 GHz 
may be useful to trace AGN activities in the central region of such 
cooling flow systems. NGC 6338 was detected at 1.4 GHz in the 
VLA FIRST survey having a flux of 57.0 mJy and radio luminosity 
of 1.37 x 10 39 erg s" 1 (Condon et al. 2002). Ratio of radio luminos- 
ity to total power of cavities is found to be ~ 10 -4 . Further, the 
entropy profile derived from the X-ray analysis is found to fall sys- 
tematically inward in a power-law fashion and is seen to flatten near 
the core. This flattening of the entropy profile at the core indicate 
toward an intermittent heating needed for maintaining its power 
when averaged over time comparable to its cooling time. Thus, all 
observational evidences in this poor cluster cooling flow galaxy in- 
dicates towards an AGN like activity working at the centre of this 
system. 
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5 CONCLUSIONS 

We have presented results based on a systematic analysis of a 
47.94 ks Chandra observation of a poor cluster cD galaxy NGC 
6338 with the major objective to study the properties of the X- 
ray cavities. The major results derived from this analysis are sum- 
marised below: 

(i) The azimuthally averaged X-ray surface brightness model as 
well as an unsharp-masked image reveals at least two , possibly 
three, X-ray cavities in the central region of NGC 6338 and are 
consistent with those seen in several other cooling flow galaxies. 

(ii) This study has detected cooler filamentary structures in X- 
ray emission spatially associating with the Ha emission line fil- 
aments and extending up to about ~7.5kpc. Dust knots seen in 
central region of optical extinction map are also found to coincide 
with these filaments and hence indicate toward a dust-aided cooling 
mechanism. 

(iii) Proximity of these cavities to the central AGN imply that 
these features are young, perhaps being inflated by jets from the 
AGN. 

(iv) Radio emission map derived from the analysis of 1.4 GHz 
VLA FIRST survey data failed to show association of radio jets 
with X-ray cavities, however, these cavities are found to be filled in 
by radio emission. High resolution, low-frequency radio observa- 
tions of this system are required to probe history of AGN outbursts 
and energy transfer. 

(v) Higher values of observed emission line ratios ([NII]/Ha , 
Ha /H/3, [OlllJ/ifg) as well as extended filamentary nature of Ha 
emission line indicate that a harder ionizing source is required to 
maintain such a high degree of ionisation. 

(vi) Ratio of radio luminosity to total power contained within 
cavities is found to be ~ 10 _4 , implying that most of the jet power 
is mechenical. 
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